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A method is presented for the identification of high-energy neutrinos from gamma ray bursts by means of a large-scale 
neutrino telescope. The procedure makes use of a time profile stacking technique of observed neutrino induced signals 
in correlation with satellite observations. By selecting a rather wide time window, a possible difference between the 
arrival times of the gamma and neutrino signals may also be identified. This might provide insight in the particle 
production processes at the source. By means of a toy model it will be demonstrated that a statistically significant 
signal can be obtained with a km 3 -scale neutrino telescope on a sample of 500 gamma ray bursts for a signal rate as 
low as 1 detectable neutrino for 3% of the bursts. 

Key words: Neutrino astronomy, gamma ray bursts, neutrino telescopes. 
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1 Introduction 

Cosmic radiation is a valuable source of informa- 
tion about various energetic astrophysical processes. 
Candidates for the production of the most energetic 
cosmic rays are Active Galactic Nuclei (AGN) and 
Gamma Ray Bursts (GRBs). 

Interactions of accelerated protons and electrons 
with the ambient photons at the acceleration site 
give rise to very energetic secondary particles. In 
particular for proton energies around the 'knee' re- 
gion of the cosmic ray spectrum, pj interactions 
yield a flux of very energetic (> 100 TeV) neutri- 
nos, mainly via production of the A resonance and 
subsequent decays. 

Various attempts^] have been made to identify 
a high-energy neutrino excess in correlation with 
satellite observations of GRBs. The performed 
searches comprise both photon-neutrino coinci- 
dence studies and investigations of so-called "rolling 
time windows". However, the former will obviously 



fail in case there exists a significant time difference 
between the arrival times of the photon and neu- 
trino fluxes, whereas the latter can only be succesful 
in case some GRBs produce multiple neutrino de- 
tections within the corresponding time windows. So 
far, no positive identifications have been reported. 

From the above it is seen that it would be prefer- 
able to use an analysis procedure that does not re- 
quire the simultaneous arrival of photons and neu- 
trinos and which also provides a high sensitivity in 
case of low signal rates. Such a method, based on 
a time profile stacking technique, is presented here 
and evaluated by means of a toy model which mim- 
ics GRB induced signals as well as (atmospheric) 
background. The only large scale neutrino telescope 
currently in operation is IceCube[^] and as such we 
use the parameters of this detector [1] as benchmark 
values for our present studies. 



* Talk presented at ICRC07, Merida, Mexico. 
1 IceCube collaboration, ICRC 2005 proceedings. 



See http://www.icecube.wisc.edu 
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2 Signal and background generation 

Our analysis procedure is based on the detec- 
tion of upgoing n tracks and as such our toy model 
only generates GRB positions homogeneously dis- 
tributed over the hemisphere opposite to the detec- 
tor location. 

For each generated burst location we define the 
satellite trigger time to be t gr b = and create a 
time window of ±1 hour around it. Observations 
with the AMANDA neutrino telescope [2] show that 
a km 3 -scale detector will observe on average 300 
upgoing muons per 24 hours due to (atmospheric) 
background, homogeneously distributed over the 
hemisphere and uniformly in time. Therefore, each 
of the above time windows will be filled with a 
number of background upgoing muon signals taken 
from a Poissonian distribution with an average 
number of 25. To conservatively account for detec- 
tor resolutions [1], Gaussian spreads of at = 10 fjs 
and a a = 1° are introduced to the arrival times and 
directions, respectively. 

Only a fraction / of the generated burst locations 
is randomly selected to yield a single upgoing fi sig- 
nal. A reasonable estimate for the possible photon- 
neutrino arrival time difference r and its spread a T 
can be obtained from the actual burst duration. 
Satellite observations [3] exhibit a mean burst du- 
ration of about 30 seconds. As such, the upgoing 
H signal arrival time of each signal burst is taken 
from a Gaussian distribution with a mean r = 30 s 
and a T = 30 s. Also for these signal muon arrival 
times and directions the corresponding detector res- 
olutions at and a a are introduced. 

In order to optimise the time bin clustering of the 
signals, the bin size should be taken to be of the or- 
der of the temporal signal spread a T . However, since 
the observed redshifts of GRBs [3] exhibit a median 
value of z = 1.9 with a spread of 1.3, cosmological 
time dilation effects have to be taken into account. 
As such, we take for the time profile bin size a con- 
servative value of 5ct t , corresponding to 150 s. 
The fraction / we keep as a free parameter in order 
to determine the sensitivity of our analysis proce- 
dure for different sizes of the GRB sample. 



3 Analysis of simulated data 

The above results in a set of identical time win- 
dows with upgoing [i arrival time recordings relative 
to the corresponding GRB trigger time. 
Stacking of all these time profiles will exhibit a uni- 
form distribution for background events. However, 
in case the data contain upgoing fi signals correlated 
with the GRBs, a clustering of data bins is expected, 
allowing the identification of correlated signals. Due 
to the cumulative character of the procedure, large 
statistics can be obtained resulting in a good sensi- 
tivity even in case of low signal rates. 

For a first investigation of the performance of 
the procedure we generated 100 GRBs in one hemi- 
sphere. All parameters were set to the values men- 
tioned above and for the fraction / we used a value 
of lO^P"! The resulting stacked time profile is shown 
in Fig. 1. 
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Fig. 1. Stacked time profile for 100 GRBs with / = 0.1. 
Further details can be found in the text. 

The corresponding stacked time profile from only 
the background signals is shown in Fig. 2. 

Comparison of the number of entries from Fig. 1 
and Fig. 2 shows that 8 of our generated GRBs in- 
duced a signal in the stacked time window. How- 
ever, due to the presence of a large background we 
are not able to identify the GRB signals on the basis 
of our observations of Fig. 1 alone. 

Restricting ourselves to an angular region of 5° 
around the GRB location will reduce significantly 
the background while preserving basically all signal 
muons. The corresponding stacked time profiles of 
our previous generation, but now restricted to an 



3 F. Halzen, et al, Astrophys. J. 527 (1999) L93. 
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Fig. 2. Stacked time profile corresponding to the back- 
ground data of Fig. 1. 

angular region of 5° around the burst location, are 
shown in Figs. 3 and 4. 
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Fig. 3. Stacked time profile for 100 GRBs with / = 0.1 
and restricted to an angular region of 5° around the 
actual burst location. 
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Fig. 4. Stacked time profile corresponding to the back- 
ground data of Fig. 3. 

Comparison of Fig. 3 and Fig. 4 allows the iden- 
tification of the GRB signals in the central bin. In 
the analysis of experimental data, however, we don't 
have access to the actual corresponding background 
distribution. As such, we need to quantify our de- 



gree of (dis)belief in a background observation solely 
based on the actually recorded signals like in Fig. 3. 

4 Bayesian assessment of the significance 

Consider a hypothesis H and some unspecified 
alternative H* in the light of some observed data D 
and prior information I. We introduce the notation 
p(H\DI) to represent the probability that H is true 
under the condition that both D and I are true. 
From the theorem of Bayes we immediately obtain 



p(H\DI) = p{H\I) p(D\HI) 
p(H*\DI) p(H*\I) p(D\HJ) 



(1) 



Introducing an intuitive decibel scale, we can ex- 
press the evidence e{H\DI) for H relative to any 
alternative based on the data D and prior informa- 
tion I as : 



e{H\DI) = 10 log 



p{H\DI) 
p(H,\DI) 



(2) 



To quantify the degree to which the data 
support a certain hypothesis, we introduce the 
Bayesian observables ip = —10logp(D\HI) and 
V>* = — 10 logp(D\H*I). Since the value of a prob- 
ability always lies between and 1, we have ip ^ 
and = 0. Together with eq. (2) we obtain 

e{H*\DI) = e(H*\I)+iP-i>* ^ e(#*|/)+V> . (3) 

In other words : there is no alternative to a certain 
hypothesis H which can be supported by the data 
D by more than ip decibel, relative to H. 
So, the value ip = —10logp(D\HI) provides the 
reference to quantify our degree of belief in H. 

In our evaluation of the stacked time profile the 
main question is to which degree we believe our 
observed distribution to be inconsistent with re- 
spect to a uniform background. This question can 
be answered unambiguously if we are able to deter- 
mine the tjj value corresponding to the uniform back- 
ground hypothesis H based on our observed stacked 
time profile. 

The process of recording uniform background sig- 
nals is identical to performing an experiment with 
m different possible outcomes {A\, ...,A m } at each 
trial. Obviously, m is in our case just the number of 
bins in the time profile, the number of trials n is the 
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number of entries and all the probabilities p k cor- 
responding to the various outcomes A k on succes- 
sive trials are equal to m _1 . As such, the probability 
p{n\ . . . n m \HI) of observing n k occurrences of each 
outcome A\. after n trials is given by the multino- 
mial distribution. This immediately yields the fol- 
lowing expression for the corresponding ip value 



Following the procedure outlined above this leads 
to the discovery sensitivities as shown in Fig. 5. 
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• (4) 



5 Discovery potential 



Evaluation of the expression of eq. (4) for the data 
displayed in Figs. 1 and 2 yields tp = 713.38 dB and 
tpbkg = 709.43 dB, respectively. Consequently, it is 
required to determine the tjj value of the correspond- 
ing background before the statistical significance of 
an observed time profile can be evaluated. 

One way to investigate background signals is to 
record data as outlined above, but with fictative 
GRB trigger times not coinciding with the actual 
t gr j). In order to have similar detector conditions for 
both the signal and background studies, background 
data may be recorded in a time span covering 1 day 
before and 1 day after the GRB observation. This 
provides at least 25 different background time pro- 
files per burst, which in turn yield the corresponding 
different stacked background time profiles allowing 
the determination of an average value 4>bkg and the 
corrresponding root mean square deviation Sbkg- 

In the case of the situation reflected by Fig. 1 
this yields $ bkg = 692.04 dB and s bkg = 21.19 dB, 
which is seen to be in excellent agreement with the 
actual background value corresponding to Fig. 2. 
Comparison of the actually observed ip value of 
713.38 dB with the reconstructed background val- 
ues immediately shows that no significant signal is 
observed. 

However, evaluation of the data corresponding to 
Fig. 3 yields tp = 218.78 dB with background val- 
ues ipbkg = 99.62 dB and Sbkg = 23.98 dB. Here a 
statistically significant signal is obtained. 

Variation of the number of GRBs allows a deter- 
mination of the minimal value of the fraction / for 
which a statistically significant signal can be ob- 
tained. Common practice is to claim a discovery in 
the case a significance in excess of 5a is obtained. 
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Fig. 5. Sensitivities corresponding to a 5a signal signifi- 
cance. 



6 Summary 

The method introduced in this report allows iden- 
tification of high-energy neutrinos from gamma ray 
bursts with large scale neutrino telescopes. The pro- 
cedure is based on a time profile stacking technique, 
which provides statistical significant results even in 
the case of low signal rates. 

The performance of the method has been inves- 
tigated by means of toy model studies based on 
realistic parameters for the future IceCube km 3 
neutrino telescope and a variety of burst samples. 
From these investigations it is seen that a 5a signif- 
icance is obtained on a sample of 500 bursts with a 
signal rate as low as 1 detectable neutrino for 3% 
of the bursts. 
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